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Abstract 

We calculate in closed analytical form the one-photon loop radiative corrections 
to muon Compton scattering ^"7 — )■ ^"7. Ultraviolet and infrared divergencies are 
both treated in dimensional regularization. Infrared finiteness of the (virtual) radiative 
corrections is achieved (in the standard way) by including soft photon radiation below 
■ an energy cut-off A. We find that the anomalous magnetic moment a/27r provides 

, only a very small portion of the full radiative corrections. Furthermore, we extend 

^ I our calculation of radiative corrections to the muon-nucleus bremsstrahlung process 

^ ■ (or virtual muon Compton scattering /x~7o — )• /x~7). These results are particularly 

relevant for analyzing the COMPASS experiment at CERN in which muon-nucleus 
bremsstrahlung serves to calibrate the Primakoff scattering of high-energy pions off a 
heavy nucleus with the aim of measuring the pion electric and magnetic polarizabilities. 
. We find agreement with an earlier calculation of these radiative corrections based on a 

Oh| different method. 

Oh' 

& PACS: 12.20.-m, 12.20.Ds, 14.70.Bh 

§s 1 Introduction and summary 
m 

^ ', At present there is much interest in a precise experimental determination of the (charged) pion 
^ ■ electric and magnetic polarizabilities, a-,,- and /3,r- Within the systematic framework of chiral 
O perturbation theory the firm prediction, — = (5. 7 ± 1.0) ■ 10~^ fm^ [1], has been obtained 
O for the (dominant) polarizability difference at two-loop order. However, this chiral prediction 
• • . is in conflict with the existing experimental determinations of — = (15.6 ± 7.8) ■ 10~^ fm^ 
from Serpukhov [2] and a,r — /^tt = (11.6 ± 3.4) ■ 10~^fm'^ from Mainz [3], which amount to 
^ values more than twice as large. Certainly, these existing experimental determinations of 
d a,r ~ /Svr raise doubts about their correctness since they violate the chiral low-energy theorem 
notably by a factor 2. 

In that contradictory situation it is promising that the ongoing COMPASS [4] experiment 
at CERN aims at measuring the pion polarizabilities with high statistics using the Primakoff 
effect. The scattering of high-energy negative pions in the Coulomb field of a heavy nucleus 
(of charge Z) gives access to cross sections for 7r~7 reactions through the equivalent photons 
method [5]. In practice, one analyzes the spectrum of bremsstrahlung photons produced 
in the reaction tt~Z — >■ 7r~Z^ in the so-called Coulomb peak. This kinematical regime is 
characterized by very small momentum transfers to the nuclear target such that virtual pion 
Compton scattering vr~7Q — )■ 7r~7 occurs as the dominant subprocess (in the one-photon 
exchange approximation). The deviations of the measured spectra (at low 7r~7 center-of- 
mass energies) from those of a point-like pion are then attributed to the pion structure 
as represented by its electric and magnetic polarizabilities [2] (taking often a,r + /^vr = 
as a constraint). It should be stressed here that the systematic treatment of virtual pion 
Compton scattering 71'% — )■ vr~7 in chiral perturbation theory yields at the same order as 
the polarizability difference an—Pn a further pion-structure effect in form of a unique pion-loop 
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correction [6, 7] (interpretable as photon scattering off the "pion-cloud around the pion"). In 
the case of real pion Compton scattering 7r~7 7r~7 this loop correction compensates partly 
the effects from the pion polarizability difference [8]. A minimal requirement for improving 
future analyses of pion- nucleus bremsstrahlung 7r~Z — )■ 7r~Z'y is therefore to include the 
loop correction predicted by chiral perturbation theory. At the required level of accuracy 
it is also necessary to include higher order electromagnetic corrections to the pion-nucleus 
bremsstrahlung process 7r~Z — )■ 7r~Z^. These radiative corrections of order a have been 
calculated recently in ref . [9] , where explicit analytical expressions for the corresponding one- 
photon loop amplitudes have been written down. Their implementation into the analysis of 
the COMPASS data is currently planned. 

The setup of the COMPASS experiment [4] gives the possibility to switch (within a short 
time) from a pion beam to a muon beam which allows for precise acceptance studies with 
a point-like (purely electromagnetically interacting) particle. In this way the muon-nucleus 
bremsstrahlung process ii~ Z ^ serves the purpose to calibrate the Primakoff scatter- 

ing events of high-energy pions off a heavy nucleus (of charge Z). Clearly, in order to preserve 
the level of accuracy one should take into account also the radiative corrections to the Bethe- 
Heitler cross section, which describes the muon-nucleus bremsstrahlung process — >■ \x~Z^ 
only at lowest order. The purpose of the present work is therefore to (re-)calculate the one- 
photon loop radiative corrections to muon-nucleus bremsstrahlung (or virtual muon Compton 
scattering). Radiative corrections to bremsstrahlung have been computed already 50 years 
ago by Fomin in refs.[10, 11]. The analytical expression given in eq.(36) of ref.[ll] (in terms of 
auxiliary functions defined by multiple integrals and derivatives thereof ) is somewhat cumber- 
some for a numerical evaluation. Also, the limiting cases discussed explicitly in refs.[10, 11] 
do not apply to the kinematics of the COMPASS experiment, in which the muon beam and 
photon energy are very large, whereas the momentum transfer to the nucleus is very small. 
In this situation it is helpful to have an independent derivation of these radiative corrections 
based on a different calculational method as employed in ref. [11]. We are using here the di- 
mensional regularization of ultraviolet and infrared divergent loop diagrams. In particular, we 
specify as a function of the five independent kinematical variables, the pertinent expressions 
for the interference terms between one-photon loop and tree diagrams summed over the muon 
spin and photon polarization states. In this detailed and explicit form our analytical results 
for the radiative corrections to virtual muon Compton scattering can be readily utilized for 
analyzing the muon Primakoff data of the COMPASS experiment. 

The present paper is organized as follows. Section 2 is devoted to the (simpler) calculation 
of the one-photon loop radiative corrections to real muon Compton scattering ^"7 — )> 
The analogous case of electron Compton scattering has been considered long ago by Brown 
and Feynman in ref. [12] and the final result for the radiative corrections has been docu- 
mented in the textbook on quantum electrodynamics by Akhiezer and Berestetskii [13]. As 
a new element we use here dimensional regularization to treat both ultraviolet and infrared 
divergencies and we write down in closed analytical form the expressions for the pertinent 
interference terms between one-photon loop and tree diagrams summed over the muon spin 
and photon polarization states. While the ultraviolet divergencies drop out in the renormal- 
izable spinor electrodynamics at work, the infrared divergencies get removed (in the standard 
way) by including soft photon radiation below an energy cut-off A. We confirm the absolute 
correctness of the final formula for these radiative corrections written in eq.(52.12) of the book 
by Akhiezer and Berestetskii [13]. In the closer analysis, we find that the muon anomalous 
magnetic moment k = ajl'n provides only a very small part of the (virtual) radiative correc- 
tions. Moreover, we compare these radiative corrections with those to (spin-0) pion Compton 
scattering and observe strong similarities concerning the energy and angular dependences. In 
section 3, we extend our calculation of radiative corrections to virtual muon Compton scatter- 
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ing /x~7o where 7q denotes the virtual (Coulomb) photon that couples to the charge 

of a heavy nucleus in the bremsstrahlung process fj^^Z — )■ fJ.'Z'j. We derive again analytical 
expressions for the pertinent interference terms between one-photon loop and tree diagrams 
summed over the muon spin and photon polarization states. Since these expressions depend 
now on five independent kinematical variables, they become rather lengthy for some terms. 
The radiative correction factor arising from the soft photon radiation off the muon (which 
cancels the infrared divergencies due to the photon-loops) is evaluated together with all its 
finite pieces in the laboratory frame. For this (soft photon) part the agreement with Fomin's 
earlier calculation [1 1] is obvious. For the (virtual) radiative corrections arising from photon 
loops the agreement between our and Fomin's calculation [11] can be verified numerically 
with good precision. 

The presentation of numerical results for the radiative corrections to the muon bremsstrah- 
lung process ijl~Z — > jJi~Z^ in the COMPASS experiment is not pursued here, since that 
requires specification and implementation of the actual experimental conditions (such as muon 
beam energy, constraints on detectable energy and angular ranges, etc.). The analytical 
results presented in this work form the basis for such a detailed study in the future. 

2 Radiative corrections to muon Compton scattering 

We start out with calculating the radiative corrections to (real) muon Compton scattering. 
The in- and out-going four-momenta of the reaction, fi~{pi) + "y{ki, ei) — >■ /i~(p2) + 7(^2, ^2) 
give rise to the (Lorentz-invariant) Mandelstam variables: 

s = (pi + h)^ = (p2 + = m\x + 1) , 

u^ipi- k^f = (P2 - kif = m\y + 1) , (1) 
^ = (Pi - P2f = {ki - k2f = -m'^{x y) , 

which obey the constraint s -\- t -\- u — 2w?, with m = 105.658 MeV the muon mass. In 
the following it will be most advantageous to work (exclusively) with the two independent 
dimensionless variables x and y defined in eq.(l) in terms of s and u. This pair of variables is 
preferred due to the presence of poles in the tree level Compton diagrams (see Fig. 1) which 
go simply as 1/x and 1/y in these variables. In its general form, the T- matrix for Compton 
scattering off a spin-1/2 particle is parameterized by six independent (complex- valued) invari- 
ant amplitudes [14], which are either even or odd under the crossing transformation s ^ u. 
Observables, like the cross section, arc then given by a quadratic form in these six amplitudes. 
We avoid the technical complication of decomposing the T-matrix of each individual diagram 
in the appropriate basis of Dirac-operators. Instead, we perform for all relevant contributions 
to the squared T-matrix directly the sums over the muon spin and photon polarization states 
(via a Dirac-trace and contractions with the metric tensor g^^). In particular, for photon- loop 
diagrams the (rf-dimensional) integration over the loop-momenta is carried out after the sum- 
mation over the spin and polarization states. With that approach in mind, the unpolarized 
differential cross section for muon Compton scattering in the center-of-mass frame, including 
the radiative corrections of relative order a, can be represented in following compact form: 




with a = 1/137.036 the fine-structure constant. Here, A and B denote the direct and crossed 
tree diagram (shown in Fig. 1) and I, II, III, IV stand for the four classes of contributing one- 
photon loop diagrams (shown in Figs. 2,3). The product symbol (8) designates the interference 
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Figure 1: Tree diagrams for muon Compton scattering. 



term between the T-matrices from two diagrams with the sums over the muon spin and 
photon polarization states aheady carried out. This multiple sum is efficiently done via a 
Dirac-trace and Lorentz-index contractions, |tr[(|^2 + m)0^'^{^i +m)0^y], where O'^'^ denotes 
the Compton tensor. The symmetrization prescription +(x ^ y) in eq.(2) generates the 
additional contributions from the crossed one-photon loop diagrams (i.e. the diagrams in 
Figs. 2,3 with crossed external photon lines). 

The advantage of working with the dimensionless variables x and y shows up already when 
evaluating the Born terms. The two tree diagrams shown in Fig. 1 lead to the following simple 
polynomial expressions: 

A(g)A = 2x-xy + A, A(g)B = B(g)A = x + y + A, B®B = 2y-xy + A. (3) 

When inserted into eq.(2), they produce the well-known Klein-Nishina cross section. In the 
center-of-mass frame the variables x and y are the related to the total center-of-mass energy 
a/s and the cosine of the scattering angle cos ^cm in the following way: 



x 



m 



y 



X 



2(1 + x) 



(2 + X + X cos 



(4) 



In the laboratory frame, on the other hand, x and y are directly proportional to the energies 
and E'^ of the incoming and outgoing photon: 
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(5) 



with the scattering angle ^lab fixed by them. In the physical region the inequalities, x > 
and —X < y < —x/{l + x) < 0, hold. The radiative corrections to muon-pair annihilation 
fj,~^fj,~ — )■ 77 are obtained by specifying x and y in eq.(2) for that (crossed) reaction. 



2.1 Evaluation of one-photon loop diagrams 

In this section, we present analytical expressions (of order a) for the interference terms be- 
tween the one-photon loops diagrams and the tree diagrams for muon Compton scattering. 
We use the method of dimensional regularization to treat both ultraviolet and infrared diver- 
gencies (where the latter are caused by the masslessness of the photon). The method consists 
in calculating loop integrals in d spacetime dimensions and expanding the results around 
d = 4. Divergent pieces of one-loop integrals generically show up in form of the composite 
constant: 

1 1 TTl 
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Figure 2: One-photon loop diagrams, (I) and (II), for muon Compton scattering. 

containing a simple pole aX d = 4. In addition, 'Je = 0.5772 ... is the Euler-Mascheroni 
number and /x an arbitrary mass scale introduced in dimensional regularization in order to 
keep the mass dimension of the loop integrals independent of d. Ultraviolet (UV) and infrared 
(IR) divergencies are distinguished by the feature of whether the condition for convergence of 
the (i-dimensional integral is c? < 4 or > 4. We discriminate them in the notation by putting 
appropriate subscripts, i.e. ^uv and ^ir. In order to simplify all calculations, we employ the 
Feynman gauge where the photon propagator is directly proportional to the Minkowski metric 
tensor g'^'^. Let us now enumerate the analytical results as they emerge from the four classes 
of one-photon loop diagrams shown in Figs. 2,3. 

Class I. These photon-loop diagrams introduce the wavefunction renormalization factor Z2 — 
1 = a(2,^//j + C,uv — 2)/27r [16] of the muon and the pertinent interference terms with the 
(direct and crossed) tree diagrams read: 



a 
2^ 



miR + ^uv-2){2-y 



y 



(7) 



l0B = ^\{2^jR + ^uv-2) 
Ztt L 



X + y + A 



X 



The additive terms, —y and —y/2, come from evaluating the Dirac-traces and Lorentz- index 
contractions consistently in d spacetime dimensions^ and expanding the results around d = 4. 
Within dimensional regularization it is crucial to follow this prescription in order to recover the 
correct low-energy limit {ki, k2 — )■ 0) for Compton scattering, namely the non-renormalization 
of the Thomson amplitude [7]. 

Class II. This photon-loop diagram involves the off-shell fermion selfenergy subtracted by the 
mass shift. The corresponding interference terms with the tree diagrams A and B read: 
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^The pole term in the ultraviolet divergence ^uv — l/{d — A) + . . . counts at this point. 
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Figure 3: One-photon loop diagrams, (III) and (IV), for muon Compton scattering. 

Class III. These photon-loop diagrams involve the half off-shell vertex correction. Both of 
them contribute equally and the pertinent interference terms with the tree diagrams read: 



III® A = -^h^uvixy - 2x - 4) + 2x^ + 4:xy + 7x + 15^ — 
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(12) 



where Li2(x) = Z^^^i n~'^x"' = x d^[^{C, — x)]~^ In^ denotes the conventional dilogarithmic 
function. A particular part of the vertex corrections (on the mass-shell) is the muon anomalous 
magnetic moment k = a/27r. Taken by its own the anomalous magnetic moment leads (via 
an appropriate derivative-coupling vertex) to the contributions: 
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(13) 



By comparing with the complete (off-shell) vertex corrections written in eqs. (11,12) one no- 
tices that there is no obvious way to localize therein the anomalous magnetic moment contri- 
butions. On-shell and off-shell effects from the photon-loop appear to be inseparably mixed 
in the interference terms III ® A and III ® B. 

Class IV. The evaluation the photon-loop box diagram is most challenging and one finds for 
its interference terms with the tree diagrams A and B: 
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+ (xy + - 2x - 12 - ^(2 + y)^ In(-x) + (^4 - 2y 
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Here, we have introduced the frequently occurring logarithmic loop function: 

L{x). 

and the auxiliary variables: 
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It is remarkable that only one single loop integral involving four propagators contributes at 
the end to the interference terms IV 0^ and IV(X>-B. This (scalar) loop integral is infrared 
divergent and it can actually be solved in terms of logarithms (L(a; + y), ln(—a;), In^ w, ln2(l — 
w)) and dilogarithms (Li2(l — w)^ Li2(ty), \Ai{h±)). As an immediate check, one verifies that 
the ultraviolet divergent terms proportional to ^uy cancel out in the sums (I+II+III) ® A 
and (I+II+III) B. 



2.2 Infrared finiteness 

In the next step we have to consider the infrared divergent terms proportional to ^//j. Inspec- 
tion of eqs. (7,8, 14,15) reveals that these scale with the Born terms A^A/x and A®B/x given 
in eq.(3). As a consequence of that feature, the infrared divergent loop corrections multiply 
differential cross section da/dflcm at leading order by a (x -H- y crossing-symmetric) factor: 
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The (unphysical) infrared divergence ^ir gets canceled at the level of the (measurable) cross 
section by the contributions of soft photon bremsstrahlung. In its final effect, the (single) soft 
photon radiation multiplies the tree level cross section da/dftcm by a factor [7, 17]: 
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which depends on a small energy cut-off A. Working out this momentum space integral by 
the method of dimensional regularization (with d > 4) one finds that the infrared divergent 
correction factor ~ ^ir in eq.(18) gets eliminated and the following finite radiative 
correction factor remains: 
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with the abbreviation W 



X 



— 4(a; + l)(a; + |/)r(l — r). We note that the terms beyond 



those proportional to ln(m/2A) are specific for the evaluation of the soft photon correction 
factor ^soft in the /i~7 center-of-mass frame with A an infrared cut-off therein. 

Having the analytical results at hand, one can show that both the virtual and the real 
radiative corrections vanish as one approaches (at fixed scattering angle ^^cm) the Compton 
threshold, \/i ^ m or a; ^ 0. This non-renormalization of the Thomson cross section 
da/dQcm — Q;^(1 + cos^ ^cm)/2m^ serves as an important check on the method of dimen- 
sional regularization where spin and polarization sums are extended from 4 to d spacetime 
dimensions. 



2.3 Results and discussion 

Wc are now in the position to present numerical results for the radiative corrections to muon 
Compton scattering /i~7 — )■ /i~7. The complete radiative correction factor is ^j-cai written 
in eq.(20) plus the sum of all interference terms (second line in eq.(2)) divided by the Born 
terms (first line in eq.(2)). First, it is important to note that our calculation (which employs 
dimensional regularization to treat both ultraviolet and infrared divergent loop integrals) 
confirms the absolute correctness of the expression written in eq.(52.12) of ref.[13] for the ra- 
diative corrections to spin-1/2 Compton scattering. Fig. 4 shows in percent the total radiative 
correction factor for four selected center-of-mass energies ^/s = (2,3,4,5)m. The detection 
threshold for soft photons has been set to the value A = 3.8 MeV.^ One observes that the 
radiative corrections become maximal in backward directions cos^cm — —1, reaching values 
up to —3% at ^/s = 5m. In Fig. 5 we show for comparison the radiative corrections to pion 
Compton scattering [7] (for an equivalent infrared cut-off of A = 5MeV). One notices that 
the overall magnitude and angular dependence of the radiative corrections are very similar 
in both cases. Interestingly, the dip in backward directions is more pronounced in the spin-0 
case. 

Next, we discuss the role of the muon anomalous magnetic moment k — a/2'jT. As a 
radiative correction it represents the on-shcll part of the vertex corrections generated by the 
photon-loop diagrams of class III (see Fig. 3). The comparison made in Fig. 6 informs us 
that the muon anomalous magnetic moment k, = a/Iir provides only a very small portion of 

^For a meaningful comparison of the radiative corrections to muon and pion Compton scattering the ratio 
of infrared cut-off A to the particle mass has to be chosen equally. 
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Figure 6: EfTect of anomalous magnetic moment k — a/27r compared to the complete photon- 
loops. 



the complete radiative corrections arising from the full set of photon-loop diagrams. Most 
noticeably, the contribution of the anomalous magnetic moment vanishes in forward directions 
cos^cm — 1) whereas the complete radiative corrections grow strongly with the center-of-mass 
energy ^/s in that region. Note that we have omitted the soft photon contribution 6^.1^^ in 
the comparison shown in Fig. 6. 

The radiative corrections in forward direction cos ^cm = 1 are of particular interest, because 
for this kinematical situation the soft bremsstrahlung contribution ^^cai^ vanishes identically. 
The soft photon emission before and after the Compton scattering process interfere in this case 
destructively with each other. After taking carefully the forward limit y —x of all terms 
calculated in section 2.1, we get for the radiative corrections to muon Compton scattering in 
forward direction: 



8 + 4x - ^ ^ , , 8 - 4x - ^ , ,1 , , 

+ 1 ReLi2(x + l) + Li2(l-x) . (21) 

The dependence of this function on the center-of-mass energy = my/x + 1 is reproduced 
by the full curve in Fig. 7. Its asymptotic behavior is {a/2n) In^ x. The dashed curve in Fig. 7 
shows for comparison the same quantity for Compton scattering off a (point-like) spin-0 
particle. Using the analytical results in section 3 of ref.[7] one deduces the expression: 

6,{xr- = -[-^ - ^ -l + -^ReU,{x+l) + _^Li2(l -x)| , (22) 

for the radiative corrections to pion Compton scattering in forward direction. The asymptotic 
behavior {2a /n) Inx in the spin-0 case turns out to be weaker than for the muon. Presumably, 
this difference has its origin in the additional magnetic coupling of the photon to a (point- like) 
spin- 1/2 particle. 
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Figure 7: Radiative correction to muon and pion Compton scattering in forward direction. 



3 Radiative corrections to muon-nucleus bremsstrahlung 

In this section we extend our calculation of radiative corrections to the muon-nucleus brems- 
strahlung process ^i^ Z — )> yrZ^. In the one-photon exchange approximation a virtual 
Coulomb photon couples to the heavy nucleus of charge Z and the bremsstrahlung pro- 
cess is governed by the virtual muon Compton scattering: \i~{p\) + 7o(9) ~^ y^'^Pi) + l^i e)- 
The in- and out-going four-momenta of that subprocess give rise to the following kinematical 
variables: 



s = (pi + (]f = {P2 + k)^ = w?'{x + 1) 



,2 



qf = m\y+l), (23) 



Ei^ mei, E2 = me2, q'^ = -m?Q , go = , (24) 

where we have also noted that the extremely small recoil energy —go = g^/2Mnuci ~ of 
the nucleus can be neglected. In the following it will be most advantageous to work with the 
five independent dimensionless variables ei,e2,x,y,Q. The unpolarized (fivefold) differential 
cross section for muon-nucleus bremsstrahlung in the laboratory frame reads [15]: 



a 



Z^a'^u\p2\ 
7r2|pi||g|4 



(25) 



with uj — ko — m(ei — 62), = "m-yGj ~ 1, |g | = ni\/Q and H is the squared amplitude 
summed over muon spin and photon polarization states. This sum is carried out most con- 
veniently via a Dirac-trace and a contraction with the (negative) Minkowski metric tensor, 
— |tr[(^2 + m)0^{j)i + m)Oi]. In that procedure the virtual Coulomb photon is given the 
polarization four- vector = (1, ) which provides in the form of Lorentz-scalar products the 
muon energies Ei—mei and £'2 = m 62 as further independent variables. 

The direct and crossed tree diagrams for virtual muon Compton scattering (named again 
A and B) are shown in Fig. 8. The open square in the left lower corner should symbolize the 
heavy nucleus to which the virtual Coulomb photon 7q couples. The contribution of these 
two tree diagrams to the squared amplitude H has the form: 



:A B\ (A B 

Htree = - + - » - + - 

.X y J \x y 



(26) 
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Figure 8: Tree diagrams for virtual muon Compton scattering. 



with 



A®B 



A® A 



2el{x-2) 



6162(0; + y 



2e\{y-2) 




where the product symbol ® designates again the interference term between (the T-matrices 
from) two diagrams with the sums over muon spins and photon polarizations already carried 
out. As it is written here in eqs. (26-29) the term -fftree reproduces to the familiar Bethe-Heitler 
cross section [15] (now reexpressed in terms of the independent variables 61, 62, x, Q). We 
note as an aside that for a spin-0 particle (e.g. a pion) -fftree reduces to a much simpler form: 



The radiative corrections of order a to virtual muon Compton scattering /i 7q — )■ /i 7 
can be summarized by the contribution: 



where I, II, III, IV, V denote the five classes of one-photon loop diagrams shown in Figs. 9,10. 
The symmetrization prescription -|-(x 61 -H- —62) supplies the additional contributions 
from the crossed one-photon loop diagrams. These are obtained by interchanging the coupling 
vertex of the left and the right photon line in the diagrams of Figs. 9,10. 

3.1 Evaluation of one-photon loop diagrams 

In this section we collect the analytical expressions for the interference terms between one- 
photon loop diagrams and tree diagrams for virtual muon Compton scattering. Such a direct 
(diagrammatic) calculation of the radiative corrections to bremsstrahlung has been deferred 
in ref.[10, 11] as being too difficult and instead the so-called mass operator method [18] has 
been used. 

Class I. These photon-loop diagrams introduce the wavefunction renormalization factor Z2 — 1 
and the pertinent interference terms with the tree diagrams [A and B) read: 




(30) 




\x y J 



\®A 




A(^ A 



2 y 

+ Ci - 6162 + - 
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(32) 
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Figure 9: One-photon loop diagrams, (I) and (II), for virtual muon Compton scattering. 
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Figure 10: One-photon loop diagrams, (III), (IV) and (V), for virtual muon Compton scat- 
tering. 



I®5 



a 
Attx 



|2(2e/ij + ^uv -2)A®B + eMx + y)-ely-elx- Q{e^ - es)' + y} . (33) 
with A® A and A ® i? given in eqs. (27,28). 

Class II. This photon-loop diagram involves the off-shell selfenergy subtracted by the mass 
shift and the corresponding interference terms with the tree diagrams read: 



II OA = - — <2^uvA®A + 



Atxx 



X + 1 V 2 ' 



4e?(8 + 8a; + x^) 



+(4eie2 - y){2x + x^ - x^) - Q{8 + lOx + x^) - Ax^ - 6x^ 



ln(— x) 

2(a: + ll 



(34) 



II® 5 



a 



Attx { ' 2(x + 1) 

+2el{x -Q)- Qx{'i + Q + X + y) + x^y 
'2e\-Q 



2eie2(2Q + x + 2Qx + y) - 2e?(Q + y) 
+ 



f (3x + 2) 



62X^(0; + 2) + 



(g + y) - 61621/ + eiQ (4 + 4x - X 



2 

+6162(16 + 20x + 4x2 + x^) + ^(x^ + x^ - lOx - 8) 



ln(— x) 
(x + 1) 



(35) 



Class III. The vertex correction to the (upper) real photon emission introduces loop functions 
that depend only on the variable x (related to s = (p2 + ^)2)- The pertinent interference terms 
with the tree diagrams (A and S) read: 



III® A 



a 



4ttx 



- 2^uvA ®A-x'^- 186? + 



9Q + 1 Aej-l-Q 
2 ^ 2(x + l) 
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X 



+|(12eie2 - 4el - 2Q - 3y - 1) + 



Q{x^ -2x-2) + x^ + 6x'^ 
x\n{—x) 



+Ax + ^e\{x^ + 8a; + 6) + {2y - 8eie2)(a;2 + 3a; + 2) 

4 



+ 



Set - 8eie2 + x + 2y + -(4ef - Q) 

X 



2(a;+ 1)2 

2n 



Li2(a; + 1) - 



TT 



6 



(36) 



Attx 



x + 1 



I - 2CuvA (8) S + |(1 + Q + 2y - 2eie2) + ^(1 + + 3?/') 

4e^ + g + 2y + 2Qy 



ef(Q-y) + eie2(g + 3)--(g + l)^ 



+^(20e2 - 4eie2 - AQ + y) + e\{Q + |/) + 6162(1 - 67/ - 3g) 
a; 



+ 



X 



y, 



^(46162 + 4g - y - 20et) + -(26162 + 26^ - 2 - g - 3y) - |(7 + 3y) 
+ 2(^1 1)2 N(Q - ^) + 2eie2(3 + g) - (g + 1)') + ^{^^iv - SQ) 
+6162(5?/ - g - 17) + 1 - 6^ + ^(3g + y) + g(3g + y- Bej)^ - 2Qy 
-eliQ + y) + ei62(2g + by-6) + 6^(5 + Q) + ^{Q-Q^- 1)] ln(-a;) 



+ 



■ (26162 (g - 2 + 4y) - 262(g + 3y) + Ael - 2y - Qy - 2y' 



+4(4ei62 + 4g - y - 206?) + 26162 - 2y - 1 



x^ 



TT 



2n 



Li2(a; + 1) - — 




(37) 



Class IV. The vertex correction to the (lower) virtual photon absorption introduces loop 
functions that depend on two variables, x and Q. The pertinent interference terms with the 
tree diagrams (^4 and B) read: 



IV (8) A 



a 
4tvx 



2A^A 



Cuv + da dbblnK(a,b;x,Q) 
Jo Jo 



+ 61(6 — 5a:) + 56162a: 



-]{6Q + 5xy) + (\la I'db } 
4 Jo JO K{a,b;x,Q) 



2ei(8 + 4g - ga: + a:y) 



+26ie2a;(2 + Q + x) - 2Q'^ - Q{A + 3x + xy) - x{2x + y + xy) 

2el{Qx + Ax-8-4:Q-xy)- 26i62a:(4 + Q + x) + Q{8 + 3x + xy) 

2g2 + x{2x + 3y + xy)] b + [26162 (g - x) - 2e\{Q + y) + 2Q + xy 



+ 



X ab 



+ 



26? (4g - ga: + xy) + (26162a; - 2g - xy){Q + a:) 



ab^ 



Ael{2 + x) + 2Q + xy 



Aei{x - 2) - 46162a; + 2g + g a^b' \ } ,(38 



IV® 5 = 



a 
Attx 







1 - 2yl (g) 5 


^uv+ [ da [ 




Jo Jo 



+ ^-^{Q + x) + '-^{12-4Q-5x-5y)-f-^-^{2 + Q + x + y) 



Jo Jo K{a,b;x,Q) 



26? (2 + Q){Q + y) + 2e\{2Q + g' + 2ga: + x^) 
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+2eie2(8 + AQ + 2x - 2y - Qy + xy) - + xy - Q^{4: + 2x + y) 
-g(4 + 3,T + x^ + 2y + xy) + [q^ - xy - 2e\{2 + Q)(Q + y) 
+2eie2(2x ^2y -2Q ^ Qy - xy) ~ 2e\{^Q + + 2x + 2Qx + x^) 



b + 



Q{2 + Q + x + y) 



+Q^{6 + 2x + y) + Q{8 + 5x + x"^ + 4y + xy) 

-2e\{Q + x) - 2eie2(Q + y)] xah^ \2e\Q{Q + y) + 2eie2(4g - + xy) 
+2e^(Q + xf -Q^ - Qx{2 + x + y) - Q^(2 + 2x + y)] + [2e^(Q + x) 

-2eie2(4 + 2g + 3a; - y) - 2el{Q + y) - Q"" - Q{2 + x + y)] + Q a%'' 



2eie2(x + y - 4) - 2el{Q + y) - 2el{Q + x) + Q{Q + 2 + x + y) \ ), (39) 



X 



with the cubic polynomial K{a, b; x, Q) = b+xa{b—l)+Qa{l—a)b in two Feynman parameters 
a,b. In the physical region, a; > and Q > 0, the integrals involving l/K{a,b;x,Q) are 
singular and in fact they become complex- valued. In order to compute their real parts (which 
are only of relevance according to eq.(31)), one performs the Feynman parameter integral 
Jq db analytically and converts the occurring (complex) logarithm \n(—ax) into ln(aa;). The 
remaining integral da is free of poles and therefore unproblematic for a numerical treatment. 
For the contributions of the crossed photon-loops, where a; > is replaced by y < 0, the two- 
parameter integrals in eqs. (38,39) are real- valued and can be evaluated numerically as they 
stand. 

On the mass-shell the vertex corrections arising from the diagrams III and IV reduce to the 
anomalous magnetic moment k — a/27r. This piece alone leads to the following interference 
terms with the tree diagrams: 



K 



A 



a 
An 



!^3el - 3eie2 + + 4) + \{7x + 5y)| , 



(40) 



a ( 



n®B = - — <^eie2(4Q + a; + 2y) -ef(Q + 2y) -e^(3g + a;) 
47ra; I 

+Q(2+| + y)+g'-|(3a: + y)|. 



(41) 



Class V. The box diagram is most tedious to evaluate. A good fraction of the occurring loop 
integrals with three and four propagators can be still be solved in closed analytical form, 
while the remaining ones have to be represented as parametric integrals. Putting all pieces 
together, we find the following rather lengthy expressions for the interference terms of the 
photon-loop diagram V with the tree diagrams (^4 and B): 
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ieliQ -A + 3x + 2y)+ 4eie2(4 - 3g - 4a; - 3y) + (a; + y)(2 + a; + y) 



+Qi2Q + 6 + 3X + 2y)l UQ + x + y) L^Q + x + y)-Q L^Q) 

^ X + y^ 



+ 

X 



8ei(2 -2x-y) + 4eie2(4x -4: + 3y) + Aelx - 2Qx -{2 + x){x + y) 
^ ^a; + y + (g + X + y)2 L'^{Q + x + y) - Q{Q + x + y) L'^{Q) 



(x + y)^ 



+(2g + x + y) V4 + g L{Q) - 2(g + x + y) V4 + g + x + y L(g + x + y) 



+ 



2ei(8 - 3g - 6x - 5y) + 2eie2(8g - 8 + 7x + 7y) + 262(2: - Q) - 2{x + y) 
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-Q{2Q + 4 + 3x + y) 



-{2ei + 2eie2 - Q - x) 



x + y 
1 



+ 



^J4 + Q + x + y L{Q + x + y)- ^Ja + Q L{Q) 
2{Q + x + y) L\Q + x + y)-2Q L^Q) 

^{x + y) + {Q + x + y)sjA + Q + x + yL{Q + x + y)- Q^Ja + Q L{Q) 

ei - 62 



x + y 



ei{x + y) + 2e2{Q + x) 



{x + yy 



A{Q + x + y) L\Q + x + y)-AQ L\Q) 
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X 



x-y + 2{Q + x + y)[^jA + Q + x + y L{Q + x + y)-^A + Q L{Q) 

2(ei - 62) 



ei{x -Q) + e2{Q + x) 



2{Q + x + y){Q + x + y-2) 



\X + y)^ 

{Q + x + y) L\Q + x + y)-Q L\Q)] + (x + y){3Q + Ax + Ay) 



-6(Q + X + y)yA + Q + x + y L{Q + x + y) + (dQ^ + 10Q{x + y) 
+3(a; + yf) ^Ja + Q L{Q) + + -{Ae^e^ - ej - 3el) + Gej 

-86162 -Ael + 2Q + 2y + 2 + x + 
, Q + 1-Ael + 46162 



i(e;-4e,e. + 3el) + ^^ 



x + 1 



3e{ + 26162 + Ae\-\-Q-y 



ln(— x) 



+ 



26^(2 - x^) + 26162(2;^ - 3x - 8) + 26^(6 + 5a;) - 2Qx - 2xy 



2 

X y 



1 



Li2(x+ 1) 



TT 



+ (Q + X-26^) 



+ [ da[ db 
Jo Jo 



-+ / (ia / dbblaK(a,b;x,Q) 
A Jo Jo 
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X 
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+^ — ^(4e2 + Q + x)ab + 2el{Q - x) ab^ -AejbH-A^A 



2+Q+x+y 
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+ [46?(1 + g + a; + y) + 46i62((5 + a; + ?/ - 2) - 126^ - Q{Q + 6 + 3x + y) 



-{A + 2x){2 + X + y) 



x{2 + Q + x + y)- 

X 



L{Q + x + y) 
^A + Q + x + y 
U2{h+) - Li2(/l_) 
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26^ (a; - 4) - 26162a; + 2Q + ^ 
ln[l + a(l - a)Q\ 
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Jo 1 + all — a 



A\n{-x) L{Q + X + y) 
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db / dc 
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+Q{2Q + 3 + 2x + 2y) 
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+Qx{2y — X — x'^ — xy) 
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x-^ 



TT 
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^2 + g+^. 

a;(/i+ - L 



^jQ + x + y L{Q + x + y) + 2Li2(w) - 2Li2(l - w) + In^ w - ln2(l - w)] 
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X 
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+ 
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K{a,b;x,Q) 

Gej + 26162(2 -Q-x-y)- 2e\{\ + Q + x + y)-Q-A-2x 
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y/A + Q + x + y 



+ 
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-2el{2Q + + 3x + Qx + 2y + 2Qy + xy + y^) + Q{2 + Q + x + yf 



-x^ - 2e\{2Q + g2 + 3x + + Qy) 



X 



io 1 + a(l 



ln[l + a(l-a)g] 



a)(g + a; + I/) 



Li2(/i+) - Li2(/i-) 4 In(-x) L(g + a; + |/) 



— h- 



V4 + g + x + y 



Jo ^0 ^0 



+ 



dc 



2ab 



G'^{a,h,c]x,y,Q) 



X 



ei(6g + 2g2 + 6a; + 2ga; + 3|/ + AQy + 2a;|/ + 2y2) + 62(8 + 2g - 2ga; 



+2x — 2x^ — y — Qy — 3xy — y'^) ei{l — ac) + 62 ac — Q{Q + 2 + x + 2y) 



+3x + y + xy + y^ 



6(ei(l 



ac) + 62 ac 



-G{a,b, c;x,y,Q) 



(43) 



Here, we have introduced the auxihary variables: 

1 



/ Q+x+y 
4+Q+x+y 



g-x-y±-^(g + a; + y)(4 + g + a; + ?/) 



(44) 

and the quartic polynomial G{a, b, c; x,y,Q) — b + a{b + c — 1 — abc) x + a{l — a)b{cy + Q) 
in three Fcynman parameters a, 6, c. The logarithmic functions I/(g) and L[Q -\- x -\- y) are 
defined by eq.(16). The (solvable) integrals /q^ dc over the third Feynman parameter c occur 
in eqs. (42,43) in the following four versions: 



/■ 

Jo 



dc 



In 



K{a,b;x,Q) 



Cdc 
Jo 



G{a,b,c;x,y,Q) N{a,b;x,y) b[l + a{l - a){Q + x + y)] ' 

b 1 

G'^(a, b, c; x, y, Q) [1 + a(l - a)(g + a; + y)]K{a, b; x, Q) ' 



(45) 
(46) 



1 -a%c 

dc ■ 



G'^{a,b,c;x,y,Q) 



N{a, b; x, y)[l + a{l - a){Q + x + y)] 
b , K{a, b; x, Q) 



In 



N'^{a,b]x,y) 6[l + a(l-a)(g + a; + y)] ' 



(47) 



2bK{a, b; x, Q) ^ K{a, b; x, Q) 

N^{a,b;x,y) ^ 6[1 + a(l - a)(g + a; + y)] 

2ab , 
N'^{a,b;x,y) N{a,b; x,y)[l + a{l - a){Q + x + y)] 

introducing a new denominator polynomial N{a, b; x, y) = x{l — ab) + y{l — a)b. It is strictly 
positive for a; > and a: + y > 0. The latter inequality imposes on the squared invariant 
momentum transfer t = {pi — ^2)^, the condition —t > Qrn? = q^. As has been discussed 
in section 3 of ref.[9] this represents no restriction for the experimentally selected events in 
the Coulomb peak, where Q < 0.1 (or even less). The term with l/K{a,b;x,Q) in eq.(46) 
requires a further analytical integration db (and conversion of logarithms into those of 
absolute values) before it can be handled numerically. 

As a good check one verifies that the ultraviolet divergencies proportional to ^uv cancel 
in the sums (I+II+III+IV) (g) A and (I+II+III+IV) (g) S. 



/ dc 
Jo 



a%c^ 



G'^{a,b,c;x,y,Q) 
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3.2 Infrared finiteness 



The infrared divergent terms proportional to C,ir in eqs. (32,33,42,43) follow again the pattern 
of the Born terms, A ^ A/x and A® B/x. The differential cross section d^a/dcudfl^dflu at 
leading order gets therefore multiplied by the (infrared divergent) factor: 



dm 

"virt 



2a 

TT 



■ A + 2Q + 2x + 2y . 
1 ^==^^=^=- L(Q + X + v) 



lIR ■ 



(49) 



On the other hand, the soft photon radiation from the muon (before and after the virtual 
Compton scattering) yields the multiplicative factor: 



(^soft ^ Oi/J, 



A-d 



d'^-H r 2pi ■ p2 



m 



m 



\i\<x 



(50) 



which includes exactly the same infrared divergence (~ ^ir) but with the opposite sign. The 
remaining finite radiative correction factor reads: 



a r r 4{2 + Q + x + y) ^ 
n\[ y/A + Q + x + y 

^^^^ ^ 



, m 
m — 
2A 



In 



ei-1 



(e2 + 



e2- 1 



_ f'^^ (2 + g + x + y)[rei + (l-T)e2] ^^ re^ + {1 - r)e2 + VR \ 
Jo 2[1 + r{l - t){Q + X + y)]VR rei + {1 - r)e2 - VR j ' 

with the abbreviation R = [rei + (1 — r)e2]^ — 1 + r(r — 1){Q + x + y). We note that the 
terms beyond those proportional to ln(m/2A) are specific for the evaluation of the soft photon 
correction factor ^goft in the laboratory frame with A an infrared cut-off therein. Alternatively, 
it can also be evaluated in the muon-photon center-of-mass frame, see herefore the expression 
^rcai'' written in eq.(24) of ref. [9]. 

Finally, it is important to compare our results with Fomin's earlier calculation of the 
radiative corrections to bremsstrahlung in refs.[10, 11]. For the soft photon contribution 6l^^^\ 
written here in eq.(51), the agreement is obvious by direct comparison with the expressions 
given in eqs. (42-45) of ref. [11]. Since Fomin uses a fictitious photon mass to regularize 
infrared divergencies our quantity is to be identified with the logarithm ln(m/m^,). The 
radiative corrections due to photon loops are summarized in Fomin's work by the term (5virt = 
-^[U^^{-x,-y,Q,e,,e2) + U'%-y,-x,Q,e2,e^)]/Uo with Uo = i^tree (see eqs.(26-29)) and 
the function 2U^^{—x, — y, ei, 62) is specified in eqs. (36-41) of ref.[ll]. We have evaluated 
the term t^yirt (setting = 0) numerically at various points in the physical region x, Q > 0, 
y < 0, ei,2 > 1 and obtained agreement with our result H\oop/Htrcc (see eq.(31)) with good 
numerical precision. Another radiative correction (called — 2«l^/7r by Fomin [11]) which 
routinely is included in the bremsstrahlung process is the (leptonic) vacuum polarization to 
the virtual photon exchange. In our notation the radiative correction factor due to muonic 
vacuum polarization reads: 



Aaj 2 5 
SttIQ ~ 6 



+ (i-^)y4TgL(Q)|, 



(52) 



with L{Q) defined by cq.(16). More significant is of course the additional radiative correc- 
tion factor (5vp,e from electronic vacuum polarization which is obtained simply through the 
substitution Q Q{m/meY = 42753(5 in eq.(52). 
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The analytical results presented in section 3 of this paper form the basis for a detailed study 
of the radiative corrections to muon- nucleus bremsstrahlung ii~Z ^ A*~-Z^7- As mentioned 
in the introduction muon-nucleus bremsstrahlung serves in the COMPASS experiment [4] at 
CERN as a calibration process for pion- nucleus bremsstrahlung 7r~Z — )■ ti~ Z'-j with the aim 
of measuring with high statistics the pion electric and magnetic polarizabilities. 

Before closing the paper let us add a few remarks. Naturally, one expects that the radiative 
corrections are approximately equal for electromagnetic processes with muons and pions. 
This expectation is to some degree confirmed by our comparison of the radiative corrections 
to muon and pion Compton scattering (of real photons) in section 2.3 (sec Figs. 4,5). The 
dominance of the real radiative corrections (due to soft photon radiation) over the virtual 
radiative corrections (due to photon-loops) for small enough infrared cut-offs A gives an 
argument in the same direction. The contribution scaling with the logarithm ln(m/2A) is 
in fact universal. 
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